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Abstract 

We introduce a new family of A-dimensional quantum superinte- 
grable model consisting of double singular oscillators of type (n, N—n). 
The special cases (2, 2) and (4,4) were previously identified as the du¬ 
als of 3- and 5-dimensional deformed Kepler-Coulomb systems with 
u(l) and su{2) monopoles respectively. The models are multisepa- 
rable and their wave functions are obtained in (n, N — n) double- 
hyperspherical coordinates. We obtain the integrals of motion and 
construct the finitely generated polynomial algebra that is the direct 
sum of a quadratic algebra (5(3) involving three generators, so(n), 
so{N — n) (i.e. (5(3) © so{n) © so{N — n) ). The structure constants 
of the quadratic algebra themselves involve the Casimir operators of 
the two Lie algebras so{n) and so{N — n). Moreover, we obtain the 
finite dimensional unitary representations (unirreps) of the quadratic 
algebra and present an algebraic derivation of the degenerate energy 
spectrum of the superintegrable model. 


1 Introduction 


The isotropic and anisotropic harmonic oscillators HI [21 El a El 0 E] are 
among the most well known maximally superintegrable systems with appli¬ 
cations in various areas of physics. However, the aspect of symmetry algebras 
generated by well-dehned integrals of motion is a complicated issue in quan¬ 
tum mechanics as recognized in early work by Jauch and Hill [1]. In the case 
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of the isotropic harmonic oscillator in A^-dimensional Euclidean space one 
can apply finite dimensional unitary representations (unirreps) and Gel’fand 
invariants and their eigenvalues of the Lie algebra su{n) to provide an al¬ 
gebraic derivation of the spectrum and the degeneracies [S]. The 3D case 
of isotropic harmonic oscillator was discussed using what is now called the 
Fradkin tensor and connected to su{3) generators [7]. There are however var¬ 
ious embeddings of the symmetry algebra, some with operators that do not 
commute with the Hamiltonian [1], in terms of su{n + l), su{n, 1) and sp{n). 
The anisotropic case has been for a long time the subject of research due to 
its various applications in nuclear physics. An analysis based on well-dehned 
integrals of motion and polynomial algebras of arbitrary order is known only 
in the 2D case [5]. 

It is also known for a long time that the harmonic oscillator in 2-dimensional 
Euclidean space is related to the 2D Kepler-Coulomb system via the so-called 
Levi-Civita or regularisation transformation uni. This transformation how¬ 
ever can only be extented to certain specihc dimensions. The 3- and the 5- 
dimensional Kepler-Coulomb systems are related to the harmonic oscillators 
in 4- and 8-dimensional Euclidean space respectively via the Kustaanheimo- 
Stiefel and the Hurwitz transformations mm and these results can be 
extended to curved spaces [13]. The Kustaanheimo-Stiefel transformations 
are connected to the so-called Stackel transformations and were used to clas¬ 
sify superintegrable systems in conformally flat spaces [m ESI US] • These are 
specihc Levi-Civita, Kustaanheimo-Stiefel and Hurwitz transformations. The 
above-mentioned connections between various models can be reinterpreted in 
terms of monopole interactions: the 4D harmonic oscillator has a duality re¬ 
lation with 3-dimensional Kepler system with an Abelian u(l) monopole ( 
also refered to as MICZ-Kepler systems) [171 [T31 ESI EEl EB EB E3] and the 
8D harmonic oscillator is in fact dual to 5D Kepler-Coulomb system with a 
non-Abelian su{2) monopole (Yang -Coulomb monopole) PH EH EH E7[ EH 

EHEDlEIlEa. 

It has been discovered that some of these duality relations can be extended 
to deformed MICZ-Kepler systems and 4D singular oscillators that are sums 
of two 2D singular oscillators. It was proven that the superintegrability and 
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multiseparability properties are preserved [331 EH EH EH EH EH] • The dual of 
the 4D singular oscillator has interesting properties and its classical analog 
has period motion [37] • If has been recognized this is also the case for the 
dual of the 8D singular oscillator [ID]. Moreover, it has been shown that 
a quadratic algebra structure exists for these two models with monopole 
interactions and their duals mm- 

In this paper, we introduce a new family of A^-dimensional superintegrable 
singular oscillators with arbitrary partition (n, N — n) of the coordinates. 
This model is a generalization of the 4D and 8D systems [3HI SOl EH] obtained 
from monopole systems via the Hurwitz transformations. However, even in 
4 and 8 dimensions only the symmetric cases (2,2) and (4,4) were studied 
[341 [3D] . We show that quadratic algebra structures exist for all members of 
the family and can be used to obtain the energy spectrum. Another main 
objective of this paper is to extend the analysis presented in [4T] to more 
complicated algebraic structures. 

The paper is organized as follow. In Section 2, we present the integrals of 
motion of the new family of superintegrable systems. In Section 3, we obtain 
the quadratic algebra and present the Poisson analog and their Casimir op¬ 
erator. We also highlight the structure of higher rank quadratic algebra and 
the decomposition. In Section 4, we generalize the realizations as deformed 
oscillator algebra, construct the Fock space and obtain the hnite-dimensional 
unitary representations (unirreps). We also using this analysis provide an al¬ 
gebraic derivation of the energy spectrum. In Section 5, we use the method 
of separation of variables in double hyperspherical coordinates and compare 
with the results obtain algebraically. In the closing section 6, we present some 
discussion with few remarks on the physical and mathematical relevance of 
these algebras. 
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2 New family of superintegrable system 


Let us consider a family of A^-dimensional superintegrable Hamiltonians in¬ 
volving singular terms with any two partitions of the coordinates (n, N — n) 


H = P- + 
2 




+ 


Cl 


xf ... -f X. 


+ 


C2 


x: 


72+1 


-f ... -f X 


2 ’ 
N 


(2.1) 


where f ^ (xi, X2,.... Xw), P ^ (pi,P2, —,Pn), ^ Edi xf, Pi ^ -ihdi and 

Cl, C 2 are positive real constants. This family of systems represents the sum 
of two singular oscillators of dimensions n and N — n. It is a generalization of 
the 4D and 8D systems obtained via the Hurwitz transformation for specihc 
case (2,2) in 4D [3Sl ES] and (4,4) in 8D jlD]. In fact model fl2.ip not only 
contains cases (2,2) and (4,4), but also cases (1,3) and (1,7), (2,6), (3,5) for 
4D and 8D respectively. An algebraic derivation of the energy spectrum has 
previously been obtained only for the two symmetric cases (2,2) and (4,4). 


The system fl2.ip has the following integrals of motion 


v2 ( N N N 
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iij=i i,j=i i=i 
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_I_ _I_ _I_ _I_ nf 

Xi -f ... i- -^n+l W ... -t- X 


■'N , 


i n N , ,2 n N 

B = AHpI- E + ^ 

^ I 1=1 i=n+l I ^ I i=l i=n+l 


Cl 


+ 


C2 


Xi ... *^72+1 ' ••• ' 

which can be verified to fulhll the commutation relation 


[H,A] = 0=[H,B]. 

The hrst order integrals of motion are given by 

Jij = XiPj - XjPi, i,j = 1,2, ....,n, 

Kij = XiPj - XjPi, i,j = n + l ,...., N. 


( 2 . 2 ) 


(2.3) 


(2.4) 

(2.5) 
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The integrals of motion A and B are associated with the separation of vari¬ 
ables in double hyper-Eulerian and double hyperspherical coordinates respec¬ 
tively. Let 

( 2 - 6 ) 

i<j i<j 

J( 2 ) and K( 2 ) represent the second order Casimir operators and fulhll the 
commutation relations 


[H, Jp,] = [H, A'pi) = [.4, Jp,] = [A, A'p,] = [B, Jp|] 

= |B.ifp)l = [Jp),irp)]=0. (2,7) 


These commutation relations can be conveniently represented by the follow¬ 
ing diagrams 
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( 2 . 8 ) 
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Ji' 


(2,9) 


The hrst diagram shows that J( 2 ) and JL( 2 ) are central elements. The sec¬ 
ond and third diagrams show that J( 2 ) and K( 2 ) are the Casimir operators 
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of so{n) and so{N — n) Lie algebras realized by angular momentum 
i,j = and Kij, i,j = n+ respectively. The models are 

minimally superintegrable and the total number of algebraically indepen¬ 
dant integrals is -|- 1. The construction of the minimally superintegrable 
systems is interesting, as research has so far mainly focused on maximally 
superintegrable systems. In the next section, we construct the quadratic 
algebra, its Casimir operator and hnite dimensional unitary representations 
which give the energy spectrum of the superintegrable systems. 


3 Quadratic algebra structure 

We present in this section quadratic algebra structure Q(3) of the super¬ 
integrable Hamiltonian systems fl2.ll) . We show how the su{N) symmetry 
algebra of isotropic harmonic oscillator is broken to Q{3)(Bso{n)(Bso{N — n). 


3.1 The isotropic harmonic oscillator and su{N) 


In this subsection we review some facts about isotropic harmonic oscillator 
and su{N). The Hamiltonian fl2.ip in the limit Ci = 0 and C 2 = 0 reduces to 
the isotropic harmonic oscillator in iV-dimension is given by 


H = 


p 

~2 




2 


(3.1) 


In this limiting case, we dehne the following operators 0 by 




a- 



iuxi), 


(3.2) 


which satisfy the commutation relation [aj,a^] = —u)h6ij. Furthermore, we 
can dehne the operators 


= 


y/2ojh 


{ai,a+} 


(3.3) 
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are hermitian and satisfy the commutation relations 

\FlF^]=5uFf-6,,Fi. (3.4) 

Hence su{N) is the symmetry algebra of the isotropic harmonic oscillator. 
We present some key briefly key elements related to an algebraic derivation 
of the energy spectrum. The su{N) symmetry algebra can be embedded in 
the non-compact algebra sp{N), dehne the operators 


rp^J _ 

-^0 - 


{at, at}’ 


\/2ujh 

and satisfying the commutators 


^ 




if;, f“i = [f;, f"] = -i.^F" - Suf;^. 


(3.5) 

(3.6) 


One can construct the Casimir operator of sp{N) in the form 

where B 2 is the Casimir operator of su{N). The eigenvalues of Q 2 and B 2 
in the representations of sp{N) and su{N) are as 


N 1 

Q2 = -j{N+-), B 2 


l{l + N){N-l) 
N 


(3.8) 


where all representations belonging to the hf-eigenstates are obtained for 
I = 1,2,..., and the energy spectrum 


E — uh{l + —), / —0,1,2,... 


(3.9) 


3.2 Quadratic Poisson algebra 


We study in this subsection the Poisson algebra of the classical version of the 
superintegrable system fl2.ll) and its Casimir operators. The system has the 
second order integrals of motion, 

I { N N 

^ ^fPj- It a^iXjPiPj 

1 * 4=1 i,j=i 



(3.10) 
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B and J( 2 ), i^( 2 ) are given by 02.31) and 02.61) respectively. It can be verified 
they satisfy 

{H,A}p = {H,B}p = 0 = {H,J^2)}p = {H,K(^2)}p, (3.11) 

where {, }p is the usual Poisson bracket. Also we can check these following 
Poisson brackets 

{A Ji2)}p = {A, K(2)}p = {B, J(2)}p = 0 = {B, K(2)}p = {J(2), A'(2)}p- 

The above commutation relations show that J( 2 ) and iP( 2 ) are second order 
Casimir operators and central elements. We now construct a new integral of 
motion from A and B as 


{A,B}p = C. (3.12) 

The integral of motion C is cubic function of momenta. After a direct but 
involving computation relying on properties of the Poisson bracket and identi¬ 
ties, we can show that the integrals of motion generate the quadratic Poisson 
algebra QP{3), 

{A,B}p = C, (3.13) 

{A, C}p = -AAB + J^ 2 )H - K( 2 )H + 2(ci - C 2 )H, (3.14) 

{B, C}p ^ 2B^ - 2H^ + - 4w'^J(2) - 4 l2^K(2) - 8w^(ci + C 2 ). 

(3.15) 

The Casimir operator of this quadratic Poisson algebra can be shown to be 
cubic and explicitly given by 

K = + 4AB^-2[J^2)H - K^2)H+ 2{ci-C2)H]B+ 16AA^ 

—2[8cj^(ci -|- C 2 ) + 4ciJ^t/(2) + 4:Uj‘^K(2) + 2H‘^'\A. (3.16) 

By means of explicit expressions as functions of the coordinates and the 
momenta for the generators A, B, C and the central elements, the Casimir 
operator becomes 

Ai = -2J^2)H^-2K^2)H^-Aci + C2)H‘^-Aj12)-AkI)+2AJ^2)K^2) 

-4a;^(ci - C 2 ) J(2) + 4:A{ci - C2)K(2) - 4cu^(ci - 02 )^. 




The quadratic Poisson algebra and the Casimir operator are useful in deriving 
the quadratic algebra and Casimir operator: the lowest order terms of h 
in quantum case coincide with Poisson analog. The hrst order integrals of 
motion generate a Poisson algebra isomorphic to so{n) Lie algebra 

{Jij j Jkl}p ^ikJjl T ^jlJik ^ilJjk ^jkJili 

for i,j, k,l = 1,n and so{N — n) Lie algebra 

{^iji ^kl\p ^ik^^jl T ^jlKik ^ilKjk ^jkKili 

for i,j, k,l = n + 1,N — n. So the full symmetry algebra is a direct sum 
of the quadratic Poisson algebra QP{3), so{n) and so{N — n) Lie algebras. 

3.3 Quadratic algebra 

We now construct integral of motion C of the quantum system from 02.21) 
and 02.31) via commutator 


[A,B] = C, (3.18) 

where C represents a new integral of motion and is a cubic function of mo¬ 
menta. The cubicness of C explains the impossibility of expressing C as a 
polynomial function of other integrals of motion, which are quadratic func¬ 
tion of momenta. After an involving but direct computations that is based 
on properties of commutators and various identities, we obtain the following 
quadratic algebra (5(3) of the integrals of motion H, A and B 

[A,C] = 2h^{A,B}-h^J^2)H + h^Ki2)H-—{8ci-8c2 

fc4 

-{N -4:){N- 2n)h‘^} H + —A^(A^ - 4)S, (3.19) 

[5, C] = -2tfB^ + 2h'^H^ - imWA + Ah^uj‘^Ji2) + AhWK(^2) 

+8h^u)‘^{ci + C 2 — —n{N — n)}. (3.20) 
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It can be demonstrated, the Casimir operator is a cubic expression of the 
generators and is explicitly given in terms of the generators {A, B and C) as 


K = 


- 2h^{A, 5^} + —{16 - N{N - + 2h^ 


+-{ 8 ci - 8 c 2 -{N - 4)(iV - 2n)h^}H 


+2h^ 8a;^{ci + C 2 —— n)} + 4a;^J(2) + 4a;^it'(2) + 2H^ 


J^2)H - Ki2)H 

B - IQh^u'^A^ 

A. 


(3.21) 


Using the realization of the integrals of motion A, B, C and the central 
element as differential operators, we can represent the Casimir operator 03.211) 
only in terms of the central elements H, J( 2 ) and 


K = 2hV(2)//" + 2h2K(2)//" + —[16 ci + 16c2-{4(A^-4) 


-{N - 2nf}h^ 


^ ^(2) - 2ra;V(2)i^(2) 


+4h"^a;"^{ci — C 2 — -^{N — A){N — n)fr}J{2) ~ ^nuj^{ci — C 2 


1 


n{N - A)H^}Ki 2 ) + (ci - C 2 )' - - n){N - 4)h^ci 


4 

— -n{N — 4)h^C2 + -n{N — n){N — 4)^“^ 


(3.22) 


This is key step in the application of the deformed oscillator algebra approach 
which relies on both form of the Casimir operators. The hrst order integrals 
of motion, that are simply components of angular momentum, generate an 
algebra isomorphic to the so{n) Lie algebra 


[t/jj, Jkil T ^jlJik ^ilJjk ^jkJil)^i 


for z, j, /c, / = 1,..., n and so{N — n) Lie algebra 

\^iji ^kl\ ^(^^ik^jl T ^jl^ik ^il^jk 


for i,j, k,l = n + 1 ,..., N — n. So the full symmetry algebra is a direct sum 
of the quadratic algebra <5(3), so{n) and so{N — n) Lie algebras. Thus the 


10 











su{N) Lie algebra generated by the integrals of motion of the iV-dimensional 
isotropic harmonic oscillators is deformed into higher rank quadratic algebra 
Q(3) © so{n) © so{N — n) for the family of superintegrable systems fl2.ip . 
The structure constants of the quadratic algebra involve three central ele¬ 
ments, the Hamiltonian and the two Casimir operators of the Lie algebras 
approached in the decomposition. 


4 Energy spectrum and unirreps 


We now consider the realizations of the quadratic algebra f (l3.18l) - (l3.2(J|) i in 
terms of deformed oscillator algebra |12],|13] dehned by 

[K,6] = -6, 66t = <h(K + l), (4.1) 

where K is the number operator and the function <h(a;) is well behaved real 
function satisfying the boundary condition 


<h(0) = 0, $(x) >0, Vx > 0. 


(4.2) 


The <h(x) is the so-called structure function. The realization of Q{?>) is of the 
form A = A(K), B = 6(K) + + p()^)b, where A{x), b{x) and p{x) are 

functions. Similar to the quadratic algebra for 2D superintegrable systems 
m, we have 


/1(N) 

B(«) 

P(«) 


8 ci - 8 c 2 + 4J(2) - 4iL(2) + (4A^ - 8n + 2nN - N‘^)h^ 

+u)^- ^ 

3.220.ni6(K + «)(1 + K + m ){1 + 2(K + u)Y' 


where w is a constant ( in fact a representation dependent constant ) to be 
determined from the constraints on the structure function. We construct 
the structure function $(x) by using the realization, the quadratic algebra 
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f fl3.18p . 03.191) . 03.201) ) and the two forms of the Casimir operator 03.211) and 


<h(a:; u, H) 


12288h^2 


64c? + 64c? - 48h^ - 32hV(2) + 16 Jj) “ 


—32 J(2)K(2) + 16/’^(2) ~ 64h^ J(2)n + QAfi?K{2)n + 48h'^n^ 
+32h^iV + 32h^ J(2)A^ — 32h^K(^2)N — ASfi^nN + 16h^ J(2)niV 
-l&h^K(2)nN - ?,2h%^N + 8h^Ar2 _ ^ 




-16c2[4(J(2) - K( 2 )) + h^{{N - 4)(2n - N) + 4(1 - 2{x + m))^}] 
-16ci[8c 2 - 4J(2) + 4Ar(2) + h^{{N - A){N - 2n) 

+4(1 - 2{x + m))'}] + 32h2[4(J(2) + K( 2 )) + h^{2n^ + {N - 2f 
—2?t,A^}](x + m) — 32h^[4(J(2) 4“ ^{2)) 4“ h^{2(?T.^ — 2) 

—2(n + 2)N + N‘^}]{x + uY — 512K^{x + u)^ + 256h‘^{x + 
x[H^-h'^{l-2{x + u)yu^]. (4.6) 


We will show how the hnite dimensional unirreps can be obtained using an ap¬ 
propriate Fock space. In 2D, we need |n, E > such that K|n, E >= n|n, E >. 
However, in our case one needs to dehne quantum numbers associated to cer¬ 
tain subalgebra chain. We use two subalgebra chains, so{n) D so{n — 1) D 
... D so{2) and so(A^ — n) D so{N — u — 1) D ... D so(2) related chains of 
quadratic Casimir operators m and 1 ^ 2 “^ can be written as 


a 


^ ^ Jij : ^ 2, ..., 72, 

(4.7) 

Q. 

Y^Kij, a = n + 2, ...,7V. 

(4.8) 


In fact, |n, i? > means |n, i?, /at, ..., /„+ 2 , In, >• Then the eigenvalues of 
J( 2 ) and iC( 2 ) are h^/„(/„+n—2) and h^ljq_ni}N-n+N—n—2) respectively. The 
constraint 04.2p can be imposed on Fock type representation of the deformed 
oscillator algebra, \n,E > with 'R\n,E >= n\n,E >, and H in ^{x,u,H) 
replaced by E. Hence by using the eigenvalues of J( 2 ), iC( 2 ) and H, the 
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structure function becomes the following factorized form: 


$(a;) 


— 12582912/i^®a;^[x + u — -(2 + mi + m 2 )] [x + u — -(2 — mi + m 2 )] 
x[x + u- ^(2 + mi - m2)][x + u- ^(2 - mi - m2 )] 


x{x + u 


—H Hu! 

2 hu) 


){x 


H fujj 
2t}w 


( 4 . 9 ) 


where h ^ m ^ = 8ci + 4J ( 2 ) + h^(n — 2 )^ and h ^ m ^ = 8c2 + 4K ( 2 ) + h ^( A ^ —n — 2 )^. 
For unirreps to be finite dimensional , we impose the following constrains on 
the structure function : 


$(p + 1 ; M , i ?) = 0 ; < F ( 0 ; n , i ?) = 0 ; < h ( x ) > 0 , Vx > 0 , ( 4 . 10 ) 


where p is a positive integer and u is arbitrary constant . We then obtain hnite 
(p + l)-dimensional unirreps . The solution of the constraints 04 . 101 ) gives the 
energy E and constant u. Thus we obtain the following possible structure 
functions and energy spectra , for ei = ± 1 , £2 = ± 1 , r] = 24576 h ^® a ;^. 

Set - 1 : 


u = 


—E -|- Huj 
2hu 


E = 2 huj{p -|- 1 -|- 


eimi e2m2 

4 


1 


< h ( x ) = pxh ^® a ;^[ 4 - I - 4 p — 4 a ; — (1 — ei ) mi - I -( 1 - I - e2 ) m2 ] 

X [4 - I - 4 p — 4 a ; -t- (1 - I - ei)mi — (1 — e2)m-2\ [4 -|- 4 p — 4 x 
-(1 - ei)mi - (1 - e2 ) m2 ] [A + Ap-Ax + {1 + ei)mi 
- h(l -h e2 ) m2][4 + Ap- 2 x + cimi e2m2 ]. 


Set - 2 : 


u = 


< h ( a ;) = 


E - j - Huj 
2hu 


E = 2 Hw{p - I - 1 - I - 


eimi -h e2m2 


), 


- a ;[4 + Ap + Ax- {I - ei)mi -h (1 -h e2 ) m2 ] 

X [4 - I - 4 p - I - 4 a ; - I - (1 - I - ei)mi — (1 — ^2)1x12] [4 -|- 4 p -|- 4 a ; 
-(1 - ei)mi - (1 - e2 ) m2][4 + Ap + Ax + {I + ei)mi 
- h(l e2 ) m2][4 + Ap + 2 x + Cimi -F e2m2 ]. 
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Set-3: 



eiTUi + €21712 


4 


), 


$(x) = r]{p +1 — x) [Ax — (1 — €i)mi — (1 — 62 ) 777 , 2 ] 

X [4a: (1 + ei)?77i - (1 - 62 ) 7772 ][4a: - (1 - 61)7771 (1 + 62 ) 7772 ] 

x[4a: (1 61)7771 -h (1 62 ) 7772 ] [2 + 2p + 2x + €imi + 627772 ]. 

The structure functions 4)(x) > 0 for 61 = 1, £2 = 1 and mi,m 2 > 0. 
In the limit Ci = 0 and C 2 = 0, this results coincide with A^-dimensional 
harmonic oscillator with the following relation among the quantum numbers 
/ = 2p In-u and the algebraic derivation using the su{N) and sp{N) 

Lie algebra and their Casimir operators and eigenvalues. Let us mention that 
the value of 77 do not play a role, only the sign needs to be taken into account 
for the constraint to obtain the finite- dimensional unirreps. 

5 Separation of variables 

Each member of the family of superintegrable Hamiltonian systems fl2.ll) 
is multiseparable and allows for the Schordinger equation the separation of 
variables in double hyper Eulerian and double hyperspherical coordinates. 
We can rewrite (12.11) into the sum of two singular oscillators of dimensions n 
and V — 77 as 


H = Hi + H2, 


(5.1) 


where 



C2 


(5.2) 


(5.3) 
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and the position vectors rf = xj + ... + + ... + The 

Schrodinger equation Hil}{r, VL) = E'l/jlr, fi) can also be written as 

Hi'ifj{ri,Q) = H2i>{r2,^l) = E2'ip{r2,Q), ( 5 . 4 ) 

where Ei and E 2 are the eigenvalues of Hi and H 2 respectively. The M- 
dimensional hyperspherical coordinates are given by 

Xi = rsin($_^^_i) sin(<^)_^^_2) • • •sin(<hi), 

X 2 = r sin($_^^_i) sin($_^^_ 2 ) • • • cos(<hi), 


x^f_i = r sin($^_i) cos(<h^_ 2 ), 

XM = rcos(<hAr-i), (5.5) 


where the J\f Xj’s are Cartesian coordinates in the hyperspherical coordinates, 
{$ 1 ,..., are the hyperspherical angles and r is the hyperradius. We 

can also introduce a double type of hyperspherical coordinates system by 
considering two copies and setting M = n and M = N — n respectively for 
the Hi component and the H 2 component. The Schrodinger equation of Hi 
in n-dimensional hyperspherical coordinates 


n — 1 d 

dvi Ti dvi 


' 1 '1 


'ilj{ri,n) = 0 , 


(5.6) 


where c( = u' = ^, E[ = ^ and the grand angular momentum 

operator h?{n) which satishes the recursive formula 


— A^(n) 




(9$ 


2 

n—1 


{n 


2) cot(<h„_i) 


d 

d^n-l 


A^(n — 1) 

sin^(<h„_i)’ 


(5.7) 


valid for n > 0 and A^(l) = 0. The separation of radial and angular parts of 
05.61) can be perfomed by setting '0(ri, 12) = i?i(ri)y(12„_i). Thus, we obtain 


d‘^Ri{ri) n—ldRi{ri) 
dr\ ri dri 


+ {2E'i - J' 


2c; 


^n{ln + n — 2) 


}/?i(ri) = 0, 


(5.8) 
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A^(n)9(S)„_,) = i„(;„ + n- 


(5.9) 


where 2) being the general form of the separation constant. Now the 

radial equation 05.81) can be converted, by setting u = arf, Ri{u) = u°‘f{u) 
and f{u) = e^“/i(M), to 


+ {2(hi + -In + ~ ~ (<^1 + — 0, 


-2u' 


2 4' 


(5.10) 


where 


- {\K2^n+ ^ 


1, n — 2,„ 1, n — 2 


-c\ 


-I 
2 "■ 


(5.11) 


Set 


AT 


2 4 


(5.12) 


Then the equation 05.101) represents the confluent hypergeometric equation 
and it gives the solution in terms of special function 




2T{Ni + 2{5i + lln + f)} 


"i ^2(S, + i(„ + 2) 

(~A^ij + —/n + —}; u) (5.13) 


In order the wavefunctions to be square integrable, the parameter A^i needs 
to be positive. Hence we obtain the discrete energy Ei of Hi from 05.12p as 


El — 2hu (jVi ^ 2 ) ’ 


(5.14) 


where ai = 2Si + In + The wave equation of H 2 in (N-n)-dimensional 
hyperspherical coordinates provides similar solution ' 0 Af 2 M 2 ) replacing ri, n, 
In, c'l, (^1 by r 2 , N — n, iN-n, C 2 , S 2 respectively in 05.131) and the energy E 2 , 
replacing Ni, ai by N 2 , a 2 respectively in 05.141) . Hence the energy spectrum 
of the A^-dimensional double singular oscillators 


E 


2hu} + 1 + 


Oil + C>i2 


(6.16) 
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where p = Ni + N 2 which coincides with the energy expression obtained 
algebraically. In fact, this multiseparability properties is also shared by its 
classical analog as the Hamilton-Jacobi can also be separated in the double 
type coordinates systems. 


6 Conclusion 


In this paper, we have extended the symmetric double singular oscillators 
in 4D and 8D to arbitrary dimensions with any partition (n, iV — n) of the 
coordinates. This provide a new family of minimally superintegrable sys¬ 
tems. As main results, we construct and obtain its realization as deformed 
oscillator algebra. We also construct the hnite dimensional unitary represen¬ 
tations which enable the algebraic derivation of the energy spectrum. This 
is compared with the results from the separation of variables method. More¬ 
over, the new family may also include duals of deformed higher dimensional 
Kepler-Coulomb systems involving nonAbelian monopoles [471IT8] . 


Our systems are generalisations of one of the four 2D Smorodinsky-Winternitz 
models [43]. Further generalisations are possible. For example. 


H 


C2 


Cl 


xl + 


+ a^ni 


4- . 


■ + ••• + 


Ca 


( 6 . 1 ) 


generalizes the A^-dimensional Smorodinsky-Winternitz model [IS] to the one 
with any partition (ni, n 2 ,..., nx) such that ni + n 2 + ■ ■ ■ + nx = N. This 
model would also be superintegrable but with a more complicated quadratic 
algebra and embedded structure seen in [421 EH] • 
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